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The cosmic electron and positron excesses have been explained as possible dark 
matter (DM) annihilation products. In this work we investigate the possible effects 
of such a DM annihilation scenario during the evolution history of the Universe. 
We first calculate the extragalactic 7-ray background (EGRB), which is produced 
through the final state radiation of DM annihilation to charged leptons and the 
inverse Compton scattering between electrons/positrons and the cosmic microwave 
background. The DM halo profile and the minimal halo mass, which are not yet well 
determined from the current N-body simulations, are constrained by the EGRB data 
from EGRET and Fermi telescopes. Then we discuss the impact of such leptonic 
DM models on cosmic evolution, such as the reionization and heating of intergalactic 
medium, neutral Hydrogen 21 cm signal and suppression of structure formation. We 
show that the impact on the Hydrogen 21 cm signal might show interesting signatures 
of DM annihilation, but the influence on star formation is not remarkable. Future 
observations of the 21 cm signals could be used to place new constraints on the 
properties of DM. 

PACS numbers: 95.35.+d,95.85.Pw,98.58.Ge 
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I. INTRODUCTION 



Recently, observations of the cosmic ray (CR) positron fraction by the PAMELA experi- 
ment 111 , the total electron and positron spectra by the ATIC m experiment, the PPB-BETS 

n nn n 

[3[ experiment, the HESS [4], |5| experiment and the Fermi |6[ experiment all show interesting 
excesses when compared with the expectations of the conventional astrophysical background. 
Especially the PAMELA data confirm the hint of excess of positron fraction previously dis- 
covered by HEAT [?], Q] and AMS {{J. Many theoretical models have been proposed to 
explain these phenomena in past years, among which a class of leptonic dark matter (DM) 
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11| or see the review 12|). 



annihilation models is of particular interest (e.g. 

It was soon recognized that if DM annihilation can occur near the solar location to 
contribute to the locally observed CR leptons, it should occur everywhere in the Universe as 
long as DM exists. Shortly after the proposal of using DM annihilation to account for the 



13|4l5|, the Milky Way halo 
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data, the 7-ray and radio emissions from the Galactic center 

20] , the dwarf galaxies |2l|-|23| , the galaxy clusters [24|, |25[ and the early Universe |26j-|28| were 
studied as possible constraints and/or future probes of the DM models. These studies have 
shown that the proposed DM scenarios are either excluded or at least strongly constrained if 
the DM halo properties derived from current numerical simulations are adopted. However, 
we should bear in mind that there are still large uncertainties in our understanding of 
the DM halo structure, especially for the small scales which are not yet probed by the 
current simulations due to the limited resolution. In the present work, we adopt the leptonic 
DM scenario implied from the CR data, and investigate what kind of constraint could be 



placed on the properties of dark matter halos [25|, [29| . This exercise is useful in two ways. 
On one hand with the improvement of computing power and astrophysical observational 
precision, we will achieve better understanding of the DM halo, if in the future our DM- 
derived constraints are in conflict with these simulations or observations, then it helps to 
falsify the DM model considered here. On the other hand, it also helps to illustrate the 
possible model uncertainties in explaining the CR observations. 

In this work we shall focus on the DM annihilation during the evolution history of the 
Universe, mainly after the epoch of cosmic recombination. The DM annihilation can produce 
electrons/positrons, which can scatter with the cosmic microwave background (CMB) pho- 
tons to produce 7-rays and contribute to the extragalactic 7-ray background (EGRB). The 
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final state radiation (FSR) from DM annihilation into charged leptons (generated directly 
from the external legs of the Feynman diagram) can also contribute to the EGRB. Using the 

y b y F ermi 31] , one can test the 
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281 ] . In addition, the energy 



observations on the EGRB by EGRET 3Jjj and most recent 
model configurations and constrain the model parameters 
of the final state particles from DM annihilation will be partially deposited in the baryonic 
gas , causing ionization and heating. This effect can be imprinted on the CMB anisotropy 



32 



-l37j. the optical depth of CMB photons 38l . |39|, and the 21 cm signal from the neutral 



hydrogen during the dark age 4CH42| . Furthermore, the heating also increases the Jeans 



mass, suppressing the accretion of gas into minihalos, and enhancing the photoevaporation 
of gas in minihalos 43[. These effects might affect the early star formation history of the 
Universe. 

Different from previous studies, in this work we will normalize the particle physics pa- 
rameters of DM particles (including; mass, cross section, branching ratios etc.) to the recent 
observed CR lepton data 44J] 1 , then calculate the EGRB and derive the constraints on the 
DM halo profile and the minimal halo mass using EGRET and Fermi data. Then we discuss 
the possible effects of DM annihilation on the ionization and heating of the intergalactic gas, 
the 21 cm signal and the early structure formation. Throughout this paper we will adopt 
the cosmological parameters from WMAP 5-yr results combined with the baryon acoustic 
oscillation and Type la supernova data: Q x = 0.228, = 0.726, fi& = 0.046, h = 0.705 
,a"8 = 0.812 and n s = 0.96 45]. The uncertainties of the cosmological parameters would 
affect the detailed result of the dumpiness factor of DM. However, compared with the large 
uncertainties of other aspects in the study, e.g., the concentration model and structure for- 
mation history of DM, the uncertainties of cosmological parameters would play a negligible 
role. 



II. EXTRAGALACTIC 7- RAY BACKGROUND 

The 7-ray emission produced in the leptonic DM annihilation generally include two 
components: the inverse Compton (IC) component through interactions between elec- 
trons/positrons and the cosmological radiation field such as the CMB, and the direct FSR 



1 We adopt the parameters fitted using fixed background approach given in Rcf. 



44| (Table I). 
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associated with the annihilation process itself. The 7-ray flux observed at redshift z is 



obtained by integrating the emissivity through evolution history of the universe 
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46|, 



ME z) - C ( 1 + ^^ 2 M r d ^ (l + ^) 3 [A 2 (^) + l] dN 

<t>{E,z)- ^ 2 ^ ^ dz R{zi) ^expl r{z,z,E)\, (1) 

where m x is the mass of DM particle, (av) is the velocity weighted annihilation cross 
section, p c = 3Hq/8ttG is the critical density of the Universe at present, H(z) = 
H -s/(n x + tt b )(l + z) 3 + Q A is the Hubble parameter. Here E' = E(l + z')/(l + z), and 
IE 7 = iwlic cO^Ifsr * s ^ e ^" ra y generation multiplicity at redshift z' for one pair of 
annihilating DM particles. Note that the cascade contribution from pair production of the 
high energy photons interacting with CMB photons is not included. This treatment is rea- 
sonable due to the following two arguments. On one hand for the DM models considered in 
this work (see below Tabled]), the energy fraction that goes into FSR photons, which should 
be easier to reach the pair production threshold, is much smaller than that which goes into 
electrons/positrons. On the other hand, because of the large dumpiness factor the domi- 
nant contribution to the EGRB for most of the redshift range interested here (e.g., z < 100) 
comes from the DM annihilation at relatively low redshifts, where the pair production is also 
not important. The optical depth t(z; z', E') of 7-ray photon with energy E', propagating 
from z' to z can be obtained as 

where E" = E'(l + z")/(l + z'), a(E,z) is the absorption coefficient. In this work we 
include the following interactions for photon propagation: the photo-ionization, photon- 
nuclei pair production, Compton scattering, photon-photon scattering and phot on-p hoton 
pair production. The absorption coefficients of these processes are adopted from 43]. Note 
that for 6 < z < 1200 the Universe is neutral, therefore the Compton scattering is negligible; 
otherwise when the Universe is ionized the photo-ionization is dropped. For redshift less 
than 6 we further include the pair production process between 7-ray photons and cosmic 
infrared background, using the optical depth of the "baseline" model in 4§j|. Finally A 2 (z) 
is a factor describing the dumpiness of the DM structures. The detailed calculation of A 2 (z) 
using structure formation theory is presented in the Appendix. 

The FSR photon spectrum for the e + e~ or ji + fi~ channel, when m x 3> m e , m^, is given 
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by |49| 



dN 
dx 



al + (l 



x) 



log 



FSR 



2 



x 



(3) 

7T X 

where % — e, fi, a « 1/137 is the fine structure constant, s = 4m 2 and x = E/m x . For the 
t + t~ channel, there is an internal bremsstrahlung radiation from the chain r — > tt° — > 7 in 
addition to the decay products. The total radiation is given by [5^ 



diV 
dx 



x 



-1.31 



(6.94x - 4.93x 2 - 0.51x 3 ) 



-4.53a; 



(4) 



FSR 



The IC component of 7-ray photons is given by 

dn 



dN 
dE 



ic 



den(e) / dE e — x F KN (e,E e ,E), 



(5) 



where n(e) is the photon number density of the background radiation as a function of 
energy e, and 4§- is the quasi-equilibrium electron/positron energy spectrum from the DM 



annihilations. Following 28[ we assume that the electrons/positrons lose their energies 



instantaneously after the production from DM annihilation, 



dn 
dK 



1 



dE' 



dE'J 



(6) 



b{E e ,z 

with the energy loss rate b(E e , z) w 2.67 x 10" 17 (1 + zf (E e /GeV) 2 GeV s" 1 . Finally the 
differential Klein-Nishina cross section F^(e, E e , E) in Eq.Q is given by 



Fk-n(c, E e , E) 



4 7 2 e 



2q lng 



2q){l-q) 



(rg) 2 (l-g) 

2(1 + Tq) 



with <jt the Thomson cross section, 7 the Lorentz factor of electron, V 



(7) 

4e7/m e , and 
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r( 



^^y. For g < l/4 7 2 or q > 1 we set F KN (e, £? e , E) = 0. 



The calculated EGRB intensity is plotted in Fig. [U Here we adopt two sets of DM particle 
parameters, as compiled in Table [U which are derived by a Markov Chain Monte Carlo 



(MCMC) fitting of t 
II) data respectively 
p Q ps 0.25 



he PAMELA+ATIC (Model I) and PAMELA+Fermi+HESS (Model 



441 ] . The cross sections derived in 44j are based on a local DM density 
0.27 GeV cm -3 , which seems to be smaller than the value derived in recent 



studies 
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521 ] . Also, 53[ shows that the dynamic measurements might systematically 



under-estimate the local density. A larger local density will result in a smaller cross section 
(by a factor of ~ 2) if we normalize the electron and positron fluxes to the observational 
excesses. Since the present work is more or less an order of magnitude study, we do not 
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expect the change of factor 2 will affect the qualitative conclusions in this work. But keep 
in mind that the signals discussed in the following should be recognized as upper limits of 
such scenarios. 

The Left and right panels are for two concentration models of DM halos: Bullock et al. 



(B01, 54j) and Maccio et al. (M08, 55[) respectively. The expected EGRB intensity is a 



function of the DM halo parameters 7 and M min (see the Appendix). In this plot we adopt 



7 = 1 (NFW profile, 56|), and varying M min from 10 7 M to 10 -6 M , corresponding to 
the observational lower limit of dwarf galaxies and theoretical lower limit of the cold DM 
scenario. It is shown that for B01 concentration the parameter settings with M m ; n = 10 -6 M Q 
are allowed by EGRET data but marginally excluded by Fermi data. For M08 concentration 
model, the extension of halos to very low mass range seems to over produce 7-rays and 
conflict with both sets of data. 

TABLE I: Mass, annihilation cross section and branching ratios of DM particles used in this work. 





m x 


(av) 


B e B^ B T 




(TeV) (10- 


-23 cm 3 s -l) 




Model I 


0.71 


0.97 


0.8 0.1 0.1 


Model II 


2.21 


9.10 


0.0 0.7 0.3 



We then scan the two dimensional parameter space M min — 7, and calculate the expected 
EGRB for each setting of parameters. By requring the calculated spectra not to exceed 
2cr errorbars of the observational data, we derive the constraints on these two parameters, 
as shown in Fig. [2j The parameter space to the topdeft side of each curve is excluded. 
The left panel is for B01 concentration model and the right panel is for M08 concentration 
model. In each panel we plot two groups of lines, corresponding to the EGRET and Fermi 
constraints respectively. We can see that the constraints are sensitive to the concentration 
models. The M08 concentration model were obtained by a power-law fitting in the high mass 
range and then extrapolated to the low mass range, while for B01 model the concentration 
of the low mass haloes shows an asymptotic flattening. So for M08 model the low mass 
haloes contribute more significantly to the 7-ray intensity than that of B01 model, hence 
the constraint is more stringent for M08 model when M m [ n is small. 

From the results of the fitting, we find that the configuration of standard cold DM haloes, 
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FIG. 1: EGRB from DM annihilation compared with the EGRET and Fermi data. Left panel is 
for B01 concentration, and the right panel is for M08 concentration model respectively. NFW halo 
profile is adopted. 




FIG. 2: Constraints on the DM halo parameters M m j n and 7 from EGRET and Fermi observations 
of EGRB. Left panel is for B01 concentration, and the right panel is for M08 concentration. 



e.g., the NFW profile (7 = 1, 56[) and M m i n down to ~ 10 6 M , is constrained by the 



Fermi data. A warm DM scenario with non-thermal production of DM particles seems to 
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be safely consistent with the data [57]. 



III. IMPLICATIONS 



A. Cosmic reionization and heating 



Now we turn to the effects of DM annihilation on the cosmic reionization and heating. 
The energy deposited in the intergalactic medium (IGM) from DM annihilation can be 
expressed as [33] 



cdmO) = f(z) 



2m\ n bo 



1 + zY ■ 2m 



x. ■ 



(8) 



where n^o is the present number density of baryon, and function f(z) represents the efficiency 



of energy transfer from DM to the IGM. In 



the efficiency f(z) is calculated in detail 



and parameterized with fitting functions, without including the DM dumpiness. However, 
it can not be directly used with the existence of DM structures, because the mean free path 
of photons for some energies (e.g., ~GeV) is large enough that the evolution of structures 
is necessary to be taken into account. In this work we follow the method given in 38| to 
calculate the energy deposition rate f(z). In short f(z) is defined as 381 ] 

JdEE(f)(E,z)a(E,z) 



(9) 



Mn^(i + s)W 

where <p(E,z) is given by Eq. (JT|). The denominator of the above equation represents the 
energy output of DM annihilation without the dumpiness, and the numerator is the energy 
absorbed by IGM. Note here the energy transfer from DM annihilation to IGM is mainly 
due to the absorption of photons, no matter what the annihilation final state is. For details 
of the results of f(z) please refer to 38|]. Our calculated results are slightly different from 



that given in 



due to different branching ratios of DM annihilation. 



In addition to the energy injection 



:rom DM annihilation, X-rays can also heat the neutral 
IGM up to a high temperature [59M62| . We assume that the X-ray emissivity is proportional 
to the star formation rate, which at high redshifts is proportional to the differential increase 
of baryon collapse fraction 62[] and is given by 



e x (z) = 1.09 x lCT 3 7xA 



p bfi (l + z) 
M Mpc" : 



63|(see also |64|) 
d/coii 



dz 



l + z)hE{z) [eV cnrV 1 ], 



(10) 



2 See also the constraints derived in the Galactic center 



29( and Virgo cluster [25 ( 
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where E(z) = ^Q m (l + zf + fi A , f x ~ 1.0 is a renormalization factor, /* is star formation 
efficiency and f co u is the collapse fraction. We take fx = 1 and /* = 10~ 3 . 

Given the energy injection rate, the evolution equations of the ionization fractions and 
the IGM temperature are modified through adding the following terms to the standard 
equations [33] 



r/°^iA e DM (» + e x Q) I - x H 1 , . 

V dz J 13.6eV 3(1 + / He ) + 

V dz J 24.6 eV 3(1 + / He ) H(z)(l + z) ' 1 j 
dT igm ^ 2[e DM (2:) + e x (z)] 1 + 2x H + / He (l + 2x He 



V dz ) 3 3(1 + / He ) + 1 J 

where /fie ~ 0.083 is the ratio of the number density of helium to that of hydrogen for the 
Universe, with the helium mass abundance Y^ e = 0.25. The ionization fraction of species A 
is defined as xa = — . These set of differential equations are solved using the RECFAST 

A n A+ +n A n & 



code [65]. 



In Fig. [3] we show the changes of the total ionization fraction x- lon = xu + /ho^Hc and 
the IGM temperature as functions of redshift, in the presence of DM annihilation. The 
four curves correspond to the maximum dumpiness enhancement allowed by EGRET or 
Fermi data, for the two DM particle models. Actually, for different combinations of the 
parameters 7 and M m ; n which lies on each of the curve in Fig. [2] the EGRB is of the same 
level. However, the impact on the ionization and heating histories could be slightly different. 
The parameters we chose in this study are compiled in Table [TTJ Here we adopt the B01 
concentration model and try to include as many low mass haloes as possible (i.e., the left 
end of each curve in Fig. [2]). 

TABLE II: The halo parameters used to calculate the reionization and heating processes. B01 
concentration is adopted. 

Model I (EGRET) Model II (EGRET) Model I (Fermi) Model II (Fermi) 
7 1.37 1.42 0.85 0.90 

M min (M ) 10~ 6 10" 6 10~ 6 10~ 6 

We can see that at z ~ 800 the effects of DM annihilation on the recombination process 
begin to appear, which lead to a higher plateau of the cosmic ionization fraction. For z < 50 
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FIG. 3: The evolution of the cosmic ionization fraction (upper panel) and the IGM temperature 
(lower panel), under the constraints from EGRB. 



the dumpiness enhancement becomes effective, which retards the decreasing trends of both 
the ionization fraction and IGM temperature. This shows that the DM annihilation may 
have obvious contributions to the ionization fraction of the Universe. However, it seems that 
the DM component is not able to fully ionize the Universe as suggested in 
to violate the EGRB observations. It is also shown in 



661 ]. in order not 



that the model which can fully 



ionize the Universe will violate the CMB measurements. 

This change of the recombination history will broaden the surface of last scattering, 
suppress the temperature fluctuations and enhance the polarization fluctuations of the CMB. 
A higher level of residual ionization fraction will also increase the scattering probability of 
CMB photons after decoupling, thus increasing the optical depth of CMB photons, which 
has been measured by WMAP 45[. The implications for the DM models from the above 
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Since these effects are not sensitive 



effects on the CMB were discussed recently in 
to the dumpiness factor A 2 (z), we will not discuss them in detail here. We point out that 
the mass and cross section fitted in 44|, and also used in the current study are consistent 
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with the WMAP constraints at 2a level 



B. 21 cm signal from neutral Hydrogen 



The redshifted 21 cm signal from neutral hydrogen can provide valuable information on 
the thermal state of gas in the dark ages. The differential brightness temperature relative 
to that of the CMB is given by 

1/2 r 



ST b = 16(1- x n ) (1 + 6) 



O02 



10 



0.3 



T7~ 



[mK], 



(14) 



where zh is the hydrogen ionization fraction as defined in the previous section, Tcmb(^) is 
the CMB temperature at redshift z, 5 is the density contrast of the gas which should be 
for the cosmic mean field 3 , Q m = Q x + f2&, and T s is the spin temperature of the neutral 
hydrogen. The spin temperature, which describes the distribution on the two hyperfine-split 
ground states of neutral hydrogen, is given by a weighted average of the CMB temperature 
and gas kinetic temperature, 

Tcmb(z) + {y a , c s + y c )T igm 



(15) 



1 + Va,eS + Vc 

The weighting factor depends on the coupling between the spin system with atoms, which are 
induced either by atom-atom, atom-proton and atom-electron collisions (y c ), or by resonance 
scattering of the Lya photons (y a , c s) which has a color temperature almost equal to the gas 
kinetic temperature. Here we adopt y c from 68J (see also the original references therein). 
The effective Lya efficiency is 

M = x e--<'+*>"*C" U + M. V 1 , (16 ) 



^lO^igm 



T- 

J 1. 



gm 



in which T+ = 0.068 K corresponds to the energy split of the hyperfine structures of neutral 
hydrogen, P 10 ~ 1.3 x 10 9 J Q with J a the specific intensity of Lya radiation, and A w = 
2.85 x 10~ 15 s" 1 is the Einstein coefficient of the hyperfine spontaneous transition (see 
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70|) 



3 This formula neglects the effect from peculiar velocity, the spherical average of peculiar velocity correction 
would give a factor 4/3 in front of 5 [67[, however, it does not affect our calculation since we consider the 
21 cm signal of the cosmic mean field, i.e. 8 = 0. 
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FIG. 4: Evolution of the 21 cm signal for the four different DM models, compared with the case 
without DM annihilation. In all cases, heating from X-ray and Lya photons from stars are included, 
see the text for details. The thin solid line corresponds to 5T^ = 0. 

We firstly consider the standard picture without the energy injection from DM annihi- 
lation. During the dark age, the scattering of residual free electrons with CMB photons 
will keep the gas temperature comparable with the CMB temperature at the initial stage. 
After kinetic decoupling of the baryons from radiation at z ~ 200, the gas temperature 
drops faster, i.e. T igm ~ (1 + z)~ 2 . The mean density of the Universe is still high at this 
time, and the coupling due to collisions is strong enough to make T s almost equal to T igm , 
so we would see an absorption 21 cm signal (c.f. the line labeled "No DM" in Fig. H]). 
This absorptional signal peaks at z ~ 100. As the redshift decreases, the density drops and 
the collisional coupling becomes less effective, resulting in the spin temperature approaches 
2cmb(z) and the 21 cm absorption signal becomes weak. However, as the stars begin to 
form, Lya photons are produced, and the spin temperature deviates from Tqmb( z ) and gets 
close to T igm again. Thus before the IGM is heated above the CMB temperature, we would 
see a second absorption peak. The location of the second peak depends on the model of 
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star formation, Lya and X-ray production rate 59(. We model the production of the Lya 
background by star formation following Ref. 711 ] . The evolution of the Lya background is 
plotted in Fig. [5] The contribution to Lya from star formation rises sharply after z ~ 40. 
This model is very simple but illustrates qualitatively the main feature of the 21 cm signal 
evolution. For our adopted model, the second absorption peak is at z ~ 14. The signal 
converts from absorption to emission at z ~ 10 due to X-ray heating. 

In the above we described the case of the global signal according to the mean density 
of gas. In addition, minihalos are formed during the late dark ages. Within minihalos 
the temperature of the gas is higher and the spin-kinetic coupling is also maintained by 
collisions as the density is higher. The minihalos could produce a positive 21 cm signal, but 
even taking into account the contribution from minihalos, the overall absorption signal is 



still apparent 



61 



72 



73|. 




FIG. 5: Evolution of the DM induced Lya intensity. For comparison, we also plot the Lya 
intensity from stars in the same panel. DM annihilation provides an early Lya background, but 
the contribution from stars increases dramatically after redshift z ~ 30. 

Now we turn to the cases with DM annihilation. The annihilation of DM has two effects 
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on the 21 cm signal. On one hand, it heats the IGM up to a higher temperature and causes 
a partial ionization, as shown in Fig. [3j On the other hand, it also produces an early Lya 
background before the formation of the first stars through de-excitation. The evolutions 
of 5Tb for different DM models as listed in Table [III are plotted in Fig. HJ The transition 
redshift at which the IGM temperature exceeds the CMB temperature is also the redshift 
at which the 21 cm signal changes from absorption to emission. 

Assuming that a fraction f a of the excitation energy goes into Lya photons, the evolution 
of Lya intensity is 

r C [1 - S H + /He(l ~ X Uc )} , , s U b 

47T 3(1 + / He ) f a ii (z) 

where ni, is the number density of baryon, v a = 2.47 x 10 15 Hz is the frequency of Lya 
photon. We adopt the energy conversion fraction f a = 0.5 {40! . 42]. The evolution of the 
DM annihilation induced Lyo; intensity is shown in Fig. [5] for different models. The DM 
annihilation produces a Lya flux even during the dark ages. As the Universe expands and 
the DM density falls, the Lya flux decreases due to the drop of annihilation rate. Thanks 
to the structure formation J a shows some enhancement for z < 80. Due to the expansion of 
the Universe, however, the flux which is proportional to the proper density of the photons 
still shows a trend of decrease, and the emission from the stars dominates at lower redshift. 

The 21 cm signal with leptonic DM annihilation contribution is shown in Fig. HI During 
the dark age, the absorption signal would not be as strong as the case without DM annihi- 
lation, because the IGM is heated to higher temperatures. Also the second absorption peak 
will become weaker due to the enhancement of DM structures. In some case there is even no 
2nd absorption peak. The transition redshift from absorption to emission will be different 
from the case without DM annihilation. However, because there is large uncertainty of the 
impact on 21 cm signal from stars, it is not easy to use the transition redshift to distinguish 
the model with DM annihilation from the standard one. We conclude that the magnitude 
of the first absorption peak will be more efficient to search for the DM annihilation signal 
and constrain the DM models. 
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C. Suppression on structure formation 

Heating by DM annihilation can also suppress the accretion of gas by haloes. The mini- 
mum mass for a halo to accrete gas from environment is the Jeans mass, 

3/2 



Mi 



^ (18) 



where /i is the averaged particle weight, G is the gravitational constant, ks is Boltzmann 



constant and m# is the mass of hydrogen atom (see 43]). The Jeans mass depends on the 
temperature of the gas. In the adiabatic model the temperature of the gas would decrease 
as redshift decreases. As a result the Jeans mass would be lower. However, with the heating 
of IGM by DM annihilations, the Jeans mass would instead increase at lower redshift, as 
shown in the left panel of Fig. [6j 

Without external heating, it is thought that the first stars (population III or Pop III) 
might be formed in DM halos with a virial temperature above 10 3 K (corresponding to halos 
with mass ~ 10 5 M at z ~ 30). In such halos the trace amount of H2 in the primordial 
gas could cool the accreted gas, which would be bound further and form stars. With the 
existence of DM annihilation, the IGM gas will be heated up to higher temperatures, and 
the corresponding Jeans mass will also become higher. However, as shown in Fig. [6] the 
Jeans mass will be lower than mass with virial temperature 10 3 K, which means the impact 
on the formation of the first stars may be weak. 

To account for the effect of DM annihilation on early star formation, we also follow the 

n 

prescription given by Oh & Haiman [74] to estimate the fraction of gas (/ gas ) left in a halo 
in the presence of an "entropy floor", i.e. heating of the gas by an external source. This 
fraction is determined by K = Kiqm/K , where Kigm = ^igm(-s) /n 2 ' 3 (z) is the entropy of 
the IGM, and K Q = T v i r /n 2 / 3 (r v i r ) is the entropy of a halo at the virial radius. We interpolate 
/gas as a function of K following j^] (see also 61]), adopting 7 = 1 (i.e. NFW profile). The 



right panel of Fig. [6] shows the evolution of the minimum mass of halos that can retain at 
least 20% of the gas. This mass is several times larger than the Jeans mass, however, it is 
still less than the virial mass supposed to form the first stars. The DM annihilation is shown 
not to significantly affect the early star formation. 
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FIG. 6: Jeans mass (left) and the minimum mass of a halo that can hold at least 20% of gas 
(right) for different DM models. The two thin solid lines correspond to the virial mass of halo with 
T v i r = 10 3 K (lower) and T v ; r = 10 4 K (upper) respectively. 

IV. CONCLUSIONS AND DISCUSSIONS 

In this work we investigate the leptonic DM annihilation, as indicated by the recent CR 
lepton data from PAMELA, ATIC, HESS and Fermi, during the evolution process of the 
Universe. We mainly focus on the period after the cosmic recombination. The mass of DM 
particle, annihilation cross section and branching ratios to three lepton flavors are adopted 
to best-fit the observational CR lepton data in the Galaxy. Then we calculate the expected 
EGRB flux from DM annihilation, including FSR and electron/positron IC emission. The 
DM dumpiness effect is included in the calculation according to the Sheth-Tormen (ST) 



halo mass function |75j . Two free parameters, the central density slope 7 and the minimum 
halo mass M min , are employed to describe the properties of the DM halo. Even so, the DM 
halo still suffers the uncertainties of the concentration. We adopt two kinds of concentration 



models in this work, B01 



54| and M08 |55|. 



The calculated EGRB results are compared with the observational data from EGRET 
and Fermi. By requiring the calculated EGRB not to exceed 2a errorbars of the data, we 
can get a conservative constraint on parameters 7 and M min . We find that for the B01 
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concentration model, the constraints on the model parameters are relatively weak, while 
for M08 concentration model in which the concentration of low mass halo is very large, 
the constraints are much stronger. Fermi data can give much stronger constraints than 
EGRET data. Even for B01 concentration, Fermi data can constrain the canonical cold 
DM prediction of halos, e.g. 7 = 1 and M m j n « 10 -6 M . Since the spectral shape of the 
DM contribution is obviously different from the observational one, we would expect that 
the DM component does not contribute significantly to the EGRB. Therefore taking into 
account the background contribution we can expect that the remaining possible DM can be 
much reduced. We suggest that a warm DM scenario with non-thermal production of DM 

n 

particles [57J might be consistent with the data safely. 

Then we investigate the effects of such models of DM annihilation on the reionization 
and heating of the IGM. In our model, in the absence of DM annihilation but with X-ray 
only, the ionization fraction and the IGM temperature till z ~ 6 are O(10 -3 ) and 0(100) K 
respectively 4 . If the DM annihilation is taken into account as an energy source, the cosmic 
ionization fraction will have a higher plateau (of the order O(10 -3 )) during the evolution 
history. Also the IGM temperature will be higher by several times for redshift 10 < z < 100. 
However, the effects of DM dumpiness on the enhancement of the energy deposition rate 
are not effective. In our treatment the energy deposition process within the evolution frame 
of the Universe is taken into account, which will result in a much lower energy deposition 
rate compared with the case with instantaneous energy deposition. We further find that 
it is impossible for the DM models considered here to fully ionize the Universe without 
conflicting with the EGRB data. 

Finally we discuss the influence of the DM annihilation on the neutral hydrogen 21 cm 
signal and the suppression on structure formation, which may connect with observations 
more directly than the ionization fraction and the IGM temperature. It is found that 
the global 21 cm signal for the models with DM annihilation differs somewhat from the 
"standard" scenario without DM annihilation. Due to the rise of the IGM temperature 

4 Note that the ionization fraction seems to be lower than that derived according to observations of Lya 
forest and the WMAP measurement of the electron scattering optical depth [76j. In this work we only 
adopt a simple description of the ordinary reionization precess and do not intend to give good description 
to Lya forest and WMAP data. In fact, there are still large uncertainties on the meansurements of mean 
transmittance of Lya flux, especially at high redshifts. A detailed model of the reionization history based 
on the Lya and WMAP observations is strongly model dependent and out of the purpose of this study. 
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during the dark age, the 21 cm absorption signal is much smaller than the "standard" case. 
Also, the Lya flux supplied by the DM annihilation would make the spin temperature be 
closer to the gas kinetic temperature at z < 40. But again due to the heating, the gas 
kinetic temperature could exceed CMB temperature earlier than the "standard" case. In 
some model the 2nd absorption peak with redshift between 10 and 20 in the "standard" 
scenario may not appear. Future 21 cm observations may have the potential to test this 
leptonic DM annihilation scenario and constrain the model parameters, despite the technical 
difficulty in practice. 

The heating of the IGM will also affect the Jeans mass, which characterizes the critical 
mass of DM halo that can accrete gas. We show that with the constraint from EGRB 
observations, the leptonic DM annihilation scenarios are not able to affect the structure 
formation process significantly. 



Appendix: Clumpiness enhancement factor of DM annihilation 



The density profile of a DM halo with mass M vir can be parameterized as 

Ps 



p(r) 



p s x p(x = r/r a ), 



(a.i; 



(r/r s )T(l +r/r s f^ 

here we adopt a free parameter 7 to represent the inner profile of DM distribution inside a 
halo. The scale radius r s is determined by the concentration parameter c vir of the halo 



:(2-7)' 



where 



3M V 



47rA vir (^)p x (z) 



1/3 



(A.2) 



(A.3) 



is the virial radius, with p x (z) = p x (l + z) 3 the DM density at redshift z, A vir (z) = (187r 2 + 
82y - 39y 2 )/(l + y), y 



1 and £l m (z) 



■771 1 . The concentration 
parameter as a function of halo mass M w - ir can be extracted from N-body simulations. For 
the low mass halos which are beyond the resolution of simulations, the concentration is 
obtained by extrapolation. In this paper we considered two extrapolation models. One is 
the analytical model developed in Ref. 54J (B01), and the other is the direct extrapolation 



of the fitting results from simulations under the WMAP5 cosmology 55[ (M08). For both 
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models we employ the redshift evolution Cvi r (z) = c vir (z = 0)/(l + z) |54J. The scale density 
p s is determined by normalizing the mass of the halo to M vir . 
The annihilation 7-ray intensity from a single halo is 

L(M vir ) = 4vr f p 2 (r)r 2 dr 



f^max 

^TT^sPs / f) 2 (x)x 2 dx 

Jo 

M vil A vir (z)p x (z)x 3 m3uX J p 2 {x)x 2 dx 



(A.4) 



3 (J p{x)x 2 dx) 2 

where x max = c vir (2 — 7). For a population of DM halos with comoving number density 
distribution dn(z)/dM v i r the total annihilation luminosity is 

L tot = J dM vir ^(l + z) 3 L(M viT ) = p 2 x (l + zf x A 2 (z), (A.5) 

in which one of the (1 + z) 3 factor comes from p x (z) term in Eq. (lA.4|) . and the other (1 + z) 3 
converts the comoving halo mass function to the physical one. In the above equation we 
define the dumpiness enhancement factor as 

,2/ a A vir(^) f dn(z) Jp 2 (x)x 2 dx 3 



A w = J dM - M ^tiw^ x - (A ' 6) 

This is used in Eq. ([!]). 

Finally we specify the halo mass function. The comoving number density distribution of 
DM haloes can be expressed as 



dn(z) p y 2A 2 a 2 r . 2 . „-, , 2 du ,. , 

— — = -^—\ 1 + iav 2 ) p \ exp (-au 2 2) , (A.7) 

where v = S c (z)/a(M v - n ), 5 c (z) = 1.68/ D(z) is the critical overdensity in spherical collapse 
model, D(z) is the linear growth factor 78] . A, a and p are constants. For (A, a, p) = 
(0.5, 1, 0) it is the Press-SchechterfPS) formula [231, and for (A, a, p) = (0.322, 0.707, 0.3) 
it is Sheth-Tormen (ST) formula 75]. In this work we adopt the ST mass function. a 2 (M vir ) 
is the average variance of density field 

a 2 (M vir ) = ^ J W 2 (kR M )P s (k)k 2 dk, (A.8) 

with the top-hat window function W(x) = 3(sinx — xcosx)/x 3 . Rm = (3M vir /47rp m ) 1 ^ 3 is 
a radius within which a mass M N - ir is contained with a uniform matter density field. The 
matter power spectrum Ps(k) is expressed as 

P 5 (k)=A s (k-Mpc)^T 2 (k), (A.9) 
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where A s is normalized using a 8 , transfer function T(k) is obtained from a fit of CDM model 
SQL 

T(q) = ln(1 2 + 3 ^ 34g) [1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 ] " a25 (A.10) 
with q = k/hT and V = f2 m /iexp[— f2b(l + \/2h/Q m )]. 
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